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CHAPTER SIX 

FORCED DISPERSIVE WAVES ALONG A NARROW CHANNEL 

Linear surface gravity waves propagating along a narrow channel display interesting 
phenomena. At first we consider free waves propagating along an infinite narrow channel. 
We give the solution for this problem as a superposition of wave modes and we illustrate 
concepts like the notion of cut-off frequency. Second, we consider a semi-infinite channel 
with forced waves excited by a wave maker located at one end of the channel. As in 
the previous case, the wave field generated by the wave maker can be described as a 
superposition of wave modes. As the wave maker starts exciting the fluid, a wave front 
develop and starts propagating along the channel if the excitation frequency is above the 
cut-off frequency for the flrst channel wave mode. If the excitation frequency is below 
the cut-off frequency for the first channel mode, the wave disturbance stays localized 
close to the wave maker, and for the particular case where the excitation frequency 
matches the natural frequency of a particular channel wave modes, there is resonance 
between this particular wave mode and the wave maker, and the wave amplitude at the 
wave maker grows with time. 

Effects of non-linearity and dissipation are not taken into account. In this chapter 
we obtain and illustrate through animations the free-surface displacement evolution in 
time along a semi-infinite narrow channel excited by a wave maker at one of its ends. 

1 Free Wave Propagation Along a Narrow Waveg- 
uide. 

We consider free waves propagating along an infinite channel of depth h and width 25. 
We adopt a coordinate system x, y, z, where x and z are in the horizontal plane and y is 
the vertical coordinate. The x axis is along the channel, the lateral walls are located at 
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z = ±b and the bottom is the plane y = —h. The free surface is located at y = rj{x, z, t), 
which is unknown. We assume irrotational flow and incompressible fluid such that the 
velocity field can be given as the gradient of a potential function (f)(x,y, z,t), where t is 
the time parameterization. The linearized boundary value problem for propagation of 
free waves is given by the set of equations 



\/'^(f){x, y, z, t) =0 for — oo < x < oo, —h < y < and — b < z < b, (1.1) 

^ =0 at y = -h, (1.3) 
dy ^ ' ^ ^ 

^ =0 at ^ = ±6, (1.4) 
oz 

and appropriate radiation conditions. This is an homogeneous boundary value problem 
that can be solved by the technique of separation of variables. First we assume that the 
free waves propagating along the channel are given as a superposition of plane mono- 
chromatic waves. Due to the linearity of the boundary value problem, we need only 
to solve it for a single nono-chromatic plane wave with wave frequency co. The time 
dependence is 



exp(— zcjt), 

and now we can write the potential function (/){x, y, z, t) and the free-surface displacement 
r]{x, z, t) in the form 

(l){x, y, z, t) = (f){x, y, z) exp{-iu}t), (1.5) 
r]{x, z,t) = r]{x, z) exp{—iu}t). (1.6) 

Now the boundary value problem given by equations (1.1) to (1.4) assume the form 
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V'^4>{x, y, z) =0 ioi — oo < X < oo, —h < y < and — b < z < b, 



(1.7) 




(1.10) 



(1.8) 



(1.9) 



where we eliminated the free surface displacement rj{x, z) and reduced the boundary 
value problem to a boundary value problem in one dependent variable, (/){x, y, z). Next, 
we apply the technique of separation of variables to solve the boundary value problem 
given by equations (1.7) to (1.10). We assume the potential function (/){x,y,z) given as 



where the possible values kz is determined by the boundary condition at the channel 
walls located at z = ±6, and the possible values of the constant k are discussed below. If 
we substitute the expression given by equation (1.11) into the boundary value problem 
given by equations (1.7) to (1.10), we obtain a Sturm-Liouville problem (one-dimensional 
boundary value problem with a second order differential equation) for the function H{y), 
which is given by the equations 



where = —k'^ + k'^. The constant A represents a set of eigenvalues, which are functions 
of the wave frequency w, of the gravity acceleration g and of the depth h. 

If we apply the boundary conditions given by equation (1.10) to the potential function 
(/){x,y,z), we realize that we can use either cos{kzz) or sm(kzz) in the expression for 
(/){x,y,z) given by equation (1.11), but with different set of possible values for the 




(1.11) 



Hyy + AH{y) = 0, 
u}^H{y) + gHy = aty = 0, 



(1.12) 
(1.13) 
(1.14) 



Hy = ai y 
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constant kz- The set of values for kz are determined by the boundary condition (1.10) 
and the choice between cos{kzz) and sm{kzz). If we consider the z dependence of the 
potential (/){x, y, z) given in terms of cos^kzz), the constant kz has to assume the values 



kzn = 7^ =1= ~7~ with n as a natural number. (1-15) 

If we consider the z dependence of the potential (/){x,y,z) given in terms of sin(A^^^2;), the 
constant kz has to assume the values 

TTITT 

kzm = ^ with m as a natural number. (1.16) 
The general form of the solution for the equation (1.12) is 



H{y) = A cosh(A(y + h)) + B sinh(A(y + h)), (1.17) 

but the boundary condition on the bottom given by the equation (1.14) implies that 
B = 0. The boundary condition at the free-surface {y = 0) gives the eigenvalue equation 
or dispersion relation 



u)'^ = gAtanh{Ah) (1.18) 

for the constant A. This implicit eigenvalue equation has one real solutions Aq and an 
infinite countable set of pure imaginary eigenvalues iAi, I = 1,2, ... . Associated with 
these eigenvalues we have the eigenfunctions 



cosh(Ao/i) 

= ^"'^-'^ + "» , with i = 1,2,... (1.20) 

COS(A;/l) 

The term ex-p{ikx){ex-p(—ikx)) in the equation (1.11) above for (p{x, y,z) represents a 
wave propagating to the right (left) if the constant k is real, or a right (left) evanescent 
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wave if A; is a pure imaginary number, or a combination of both if k is complex. We 
label the constant k as the wavenumber. Since, we are interested in free propagating 
waves, we need the constant k to be a real number. The value of this constant is given 
in terms of the constants A and A;^, according to the equation 



e = - kl, (1.21) 

where the possible values of A;^ are given by the equations (1-15) and (1.16). The possible 
values of A are solutions of the dispersion relation given by the equation (1.18). Since 
we want A; as a real number, this excludes the imaginary solutions of the equation (1.18), 
so we can write the equation above in the form 



L _a2_,2 
""ti — -''^o "'zrf) 



(1.22) 
(1.23) 



where we appended the indexes n and m to the constant k to make clear its dependence 
on the eigenvalues k^n and k^m- 

Now we can write the potential function (j){x, y, z) in the form 



fr . / , \ 7-, \n COSh(Ao(?/ + /l)) 

(p{x, y,z) = } < [Am exp{ikmX) + Bm exp{-ikmX)\ — sm{k^rnZ) 

^-^ I cosh(Ao/ij 



m=— oo 

oo 



+ ^ S [-4„ exp{iknx) + exp{-iknx)] 



cosh(Ao(?/ + /0) 
cosh(Ao/i) 



COs{kznZ) 



(1.24) 



and the free-surface displacement r]{x, z) is given by the equation 



ri{x,z) 



N n / ■, cosh(Ao(?/ + /i)) . . 

[Am exp{zkmX) + Bm exp{-zkmX)) — sm{k^mZ) 

cosh(Ao/ij 



E 



cosh(Ao(?/ + /i)) 
[An exp[iknX) + Bn exp(-zA;„a;)) — cos[k^nZ) 

COSil( i Vq/2. 1 



(1.25) 
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where the value of the constants A^, An, and Bn are specified by the appropriate 
radiation conditions. 

According to the value of k^m or k^n, the constants k^ and kn in the equations (1-24) 
and (1.25) may be real (propagating wave mode) or pure imaginary numbers (evanescent 
wave mode). If we fix the value of kzm or kzn (fix the value of m or n), for a given depth 
h, we can vary the wave frequency co such that Aq > kzm{kzn) or Aq < kzm{kzn)- When 
^0 > kzm{kzn), km{kn) is a real number and we have a propagating wave mode, but 
when Aq < kzm{kzn) we have that km{kn) is a pure imaginary number and the wave 
mode associated with this value of k is evanescent. So, the wave frequency value where 
kzm = ^o{kzn = ^o) is Called the cut-off frequency for the mth (nth) wave mode. 

Next, we plot the dispersion relation given by equation (1.18) as a function of the 
wavenumber k and the depth h for various values of the eigenvalues kzm (sine wave 
modes in the z coordinate) in the figures 1 and 2. As the value of kzm increases (value 
of m increases), the wave frequency assume larger values for the considered range of the 
wavenumber k. The wave frequency value at A; = for a given kzm (given m) is the cut- 
ojf frequency for the wave mode associated with the eigenvalue kzm- For a fixed value 
of kzm, frequencies below the cut- ojf frequency implies in pure imaginary wave numbers 
and the associated wave mode is exponentially decreasing (evanescent) or exponentially 
growing. Wave modes associated with pure imaginary wave numbers do not participate 
in the superposition leading to free waves solutions. According to figures 1 and 2, the 
higher the wave frequency, the higher the number of wave modes participating in the 
superposition leading to free waves solutions. 

Another way to see that the wave modes associated with imaginary wave numbers 
(wave below the wave mode cut- ojf frequency) do not propagate is through the wave 
mode group velocity. In figures 3 and 4, we plot the group velocity for the first 10 
wave modes associated with the eigenvalues kzm ("^ from to 9). For wave frequencies 
above the cut-off frequency, the considered wave mode (fixed value of kzm) has a real 
wavenumber k and non-zero group velocity, as we can see through figures 3 and 4. As the 
wave frequency approaches the cut- ojf frequency, the group velocity of the considered 
wave mode approaches zero, according to figures 3 and 4. At the cut- ojf frequency of 
the considered wave mode, its group velocity is zero and no energy is transported by 
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wavenumber (radians per meter) 

Figure 1: Wave frequency as a function of the wavenumber k for various values of the 
eigenvalue k^m and water depth h = 100 meters. 
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wavenumber (radians per meter) 

Figure 2: Wave frequency as a function of the wavenumber k for various values of the 
eigenvalue k^m and water depth h = 0.1 meters. 
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this wave mode for wave frequencies at or below the wave mode cut-off frequency. 

According to figures 3 and 4, the group velocity for each wave mode has a maximum 
value, which decays as the value of k^m increases (value of m increases). The first wave 
mode (sine wave mode with kz = 0) has the largest maximum group velocity, and since 
its cut-off frequency is zero, we can have free propagating waves for any wave frequency 
for the channel specified by its depth h, its width 2b and the gravity acceleration g. 
Above, we looked at the wave modes with sine dependence in the z coordinate. For the 
wave modes with cosine dependence in the z coordinate, the minimum absolute value 
of the eigenvalue kzn is larger than the minimum absolute value for the eigenvalues 
kzm, which is zero. Therefore, for any wave frequency we have free waves propagating 
along the channel. For the cosine wave modes there is a minimum cut-off' frequency. 
Propagation of this type of wave mode is possible only for wave frequencies above their 
minimum cut-off frequency. 

2 Forced Wave Propagation Along a Narrow Waveg- 
uide. 

Now we consider forced waves propagating along a semi-infinite channel with the same 
depth h and width 2b as the channel in the previous section. The semi-infinite channel 
has a wave maker at one edge of the channel, which generates wave disturbances that 
may or may not propagate along the channel. The solution for the forced waves is given 
as a superposition of wave modes. The same wave modes we obtained in the previous 
section. Evanescent wave modes are also part of the solution in this case. They stay lo- 
calized close to the wave maker and describe the local wave field. For a mono-chromatic 
excitation, the wave modes with cut-off frequency below the excitation frequency con- 
stitute the propagating wave field, and the wave modes with cut-off frequency above the 
excitation frequency are evanescent and stay localized close to the wave maker. Their 
superposition gives the evanescent wave field. 



10 




Figure 3: Group velocity as a function of the wavenumber k for various values of the 
eigenvalue k^m and water depth h = 100 meters. 
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Figure 4: Group velocity as a function of the wavenumber k for various values of the 
eigenvalue k^m and water depth h = 0.1 meters. 
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3 Initial Boundary Value Problem. 

We consider the same coordinate system used in the previous section. The wave maker is 
located at a: = and the channel lays at a: > 0. The linearized boundary value problem 
for the forced waves is similar to the boundary value problem for the free waves problem. 
The difference is the boundary condition describing the effect of the wave maker and 
the fact that the channel is now semi-infinite. The linear boundary value problem for 
forced waves is given by the set of equations 



V^<^(a;, y, z,t) =0 iov < x < oo, -h < y < and - b < z < b, (3.26) 
d'^cl) del) 



+ =0 aty = 0, (3.27) 
at y = -h, (3.28) 



dt'^ dy 
dy 



^=0 at z = ±b, (3.29) 

oz 

^ =^F{z)G{y)f{t) on X = 0, (3.30) 

and the free surface displacement 77(2;, z, t) is related to the potential function (j){x, y, z, t) 
according to the equation 



1 d(j) 

r]{x, z, ^) = -~"^(^' 0' ^' (3-31) 

The function f{t) is a known function of time. Actually, we chose an harmonic excita- 
tion, so we have 



f{t) = cos{u}t), (3.32) 

where co is the excitation frequency. We need also to consider initial conditions for the 
boundary value problem above. They are given by the equations 



(t){x,y,z,0) =0, 
(l)t{x,y,z,0) =0, 



(3.33) 
(3.34) 
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where the initial condition (3.34) is equivalent to have a still free surface at t = 
{rj{x,z,0) = 0). Next, we solve the initial boundary value problem, which is dis- 
cussed in the next section. 

3.1 Solution of the Initial Boundary Value Problem. 

The first step to solve the initial boundary value problem given by equations (3.26) 
to (3.30) is to apply the cosine transform in the x variable. This results in a non- 
homogeneous Helmholtz-like equation for the potential function under homogeneous 
boundary conditions. Since the resulting equation is non-homogeneous, the solution 
is given as the superposition of the solution for the homogeneous part of the problem 
plus a particular solution that handles the non-homogeneity. To solve the associated 
homogeneous problem, we use the method of separation of variables as in the previ- 
ous section. The solution of the homogeneous problem is given as a superposition of 
modes in the y and z variables. The particular solution is obtained using the homo- 
geneous solution through the method of variation of the parameters. The constants of 
the homogeneous solution are obtained by applying the boundary conditions to the full 
solution (homogeneous plus particular solutions). Next, we discuss in detail the steps 
outlined above. 

We consider the cosine transform pair 



If we apply the cosine transform (3.36) to the second partial derivative of the potential 
function (l){x,y,z,t) with respect to the x variable, we have that 




(3.35) 



and 




(3.36) 





(3.37) 
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since we assumed that — > and — > as a: — > 00. The term (f)x{0, y, z, t) is specified 
by the boundary condition at a; = and given by equation (3.30). Next, we apply 
the cosine transform to the initial boundary value problem given by equations (3.26) to 
(3.30). This results in the set of equations 



<j)yy + <j)zz - k'^4> = 4>x{0, y, z, t) = —F{z)G{y) cos(wt), (3.38) 

4>tt + g(t)y =0 on y = 0, (3.39) 

=0 on y = -h, (3.40) 

4 =0 on ^ = ±6, (3.41) 

with the initial conditions given by equations (3.33) and (3.34) written in the form 

<^(fe,y,2;,0)=0, (3.42) 

^t(A;,y,^,0)=0. (3.43) 



This is a non-homogeneous initial boundary value problem for the function (f)(k,y,z,t) 
(cosine transform of (l){x,y, z,t)). Our strategy to solve this initial boundary value 
problem is to find the general form of the solution of the homogeneous part of the 
initial boundary value problem given by equations (3.38) to (3.41) plus a particular 
solution for the non-homogeneous part of this initial boundary value problem. To find 
the value of the constants of the homogeneous part of the solution, we apply the initial 
and boundary conditions to the full solution (homogeneous plus particular). Next, we 
consider the homogeneous part of the initial boundary value problem for </), which is 
given as the superposition of wave modes obtained in the previous section. So, the 
solution of the homogeneous problem is similar to the one given by equation (1.24). 
The solution for the homogeneous problem is 

+00 

(j>H=Yl {[Mk,t)cosh{Ar,{y + h)) + Br,{k,t)smh{Ar,{y + h))]cos{k,nz)} 

n=—oo 

+00 

+ X] {[^rn{k, t) cosh{A^{y + h)) + D^{k,t)smh{ A^{y + h))]sm{k^rnz)} , 
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where = k'^ + k'^^, = + and kzn and kzm are given respectively, in equations 
(1.16) and (1.15). As we mentioned before, the general solution is given as a superpo- 
sition of the homogeneous solution (pH plus a particular solution. We assume that the 
particular solution has the form 

+ 00 

4>p=^ |Ji„(/c,y,t)cosh(A„(y + /i)) + E„(/c,y,t)sinh(A„(y + /i)) cos{kznz)^ 

n=—oo 

+ 00 

+ ^ ^^^Crnik,y,t)cosh{Arn{y + h)) + Drn{k,y,t)s'mh{Arn{y + h)) s'm{kzmz)^ . 



We substitute the potential (/)p in the non-homogeneous Helmholtz equation (3.38) in 
the y and z variables. We also impose that 



d(l)p 
dy 



+ 00 

{An i„sinh(A„(y + /i)) + ^„cosh(A„(y + /i)) cos{kznz)^ 



(3.44) 



+ Yl l^"^™ sinh(A^(y + h)) + cosh(A^(?/ + h)) 



sm{kzmz) 



The procedure above results in the set of equations for the amplitudes An,Bn,Cm and 



{An)y cosh(An(|/ + h)) + {Bn)y sinh(An(|/ + h)) =0, 
{Cm)y cosh(A„(y + h)) + {Dm)y sinh(A„(y + h)) =0, 



(3.45) 
(3.46) 



A„{(i„)2,sinh(A„(y + /i)) + (^„)^cosh(A„(y + /i))} =-^G{y) cos{u}t)F^, (3.47) 
A™ I (C'm)^, sinh(A„(y + /i)) + {Dm)y cosh(A„(y + /i)) } =^^(2/) cos(a;t)F„, (3.48) 
where 



Fm= F{z) sm{kzmz)dz, 
J-b 

Fn= F{z) cos{kznz)dz. 
J-b 



(3.49) 
(3.50) 
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If we solve the set of equations above and integrate with respect to the y variable 
from —h to 0, we obtain the following expressions for the amplitudes An,Bn,Cn and 



D„, which follows: 



Cri 



Aoj 



cos{u}t)FnGniy), 



Au 



cos{iot)FnHn{y), 



Aoj 
6^ 



COs{u)t)FmGmiy), 



Au 



■co^{u}t)FmHm{y), 



(3.51) 
(3.52) 
(3.53) 
(3.54) 



where the functions Gn{y), Hn{y),Gm{y) and Hm{y) are given by the equations 



/y 
G{p) smh{An{p + h))dp, 
■h 
ry 

Hn{y) = / G{p) cosh(A„(p + h))dp, 

J-h 

/y 
G{p) sinh(A„(p+ /i))o!p, 
-h 

/y 
G{p) cosh{A^{p + hj)dp. 
■h 

Now, the total solution (/)(A;, y, z) can written in the form 



(3.55) 
(3.56) 
(3.57) 
(3.58) 



E 



An - cos{u)t)FnGn{y) ) cosh(A„(y + h)) 



Aco 



+ \Bn + TT?- cos{u}t)FnHn{y) sinh(A„(y + h)) 



E 



Aoj 



Cm - cos{u}t)FmGm{y) ) cosh(Am(y + h)) 



+ \Dm + 



Aoj 
6^ 



cos{u}t)FmHm{y) sinh(A„(y + h)) 



(3.59) 



In the expression above we still need to obtain the constants A^., B^, and of the 
homogeneous part of the solution. To do so, we apply the boundary conditions (3.39) 
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at y = and (3.40) at y = —h. The boundary condition at y = —h, given by the 
equation (3.40), implies that Dm = 0{Bn = 0). The boundary condition at y = gives 
the equation 

A F 

{Ar,)tt + c/A„tanh(A„/z)A„ =I^T^ {u^ cos (cot) [i7„(0) tanh(A„/0 - ^^(0)] 

(3.60) 

+9AnC0 cos{cot) [Gn{0) tanh(A„/i) - i/„(0)]} . 

We also obtain a similar equation for Cm- This is a non-homogeneous second order dif- 
ferential equation in time for the amplitude An- Its solution is given as the superposition 
of the solution of the homogeneous part of the equation plus a particular solution which 
satisfies the non-homogeneous term in the equation (3.60). The homogeneous solution 
is given as 

{An{t))H = A cos(n„t) + B sm{nj) (3.61) 
with = ^A„ tanh(A„/i). We assume the particular solution given in the form 



{An{t))p = A{t)p cos{nt) + B{t)p sin(fit). (3.62) 

We impose that 



^(-4„(t))p = Qn |-i(t)psin(nt) + ^(t)pcos(nt)| . (3.63) 

If we substitute the form of the particular solution, given by equation (3.62) into the 
governing equation (3.61) and take into account the assumed form for ^{An{t))p, given 
by equation (3.63), we obtain for the amplitudes A{t)p and B{t)p the expressions 



^ lV'(w,n„,/l) JcOs[(fin - W)t] C0s[(fin + W)t] . 

^ 10(w,Q„,/i) Jsin[(fi„ - w)t] sin[(fi„ + w)t] . 
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where 



AR 



xP{oj, Qn, h) =^-^ {o;^ [i/„(0) tanh(A„/i) - ^^(O)] 
-^gKui [Gn{Q) tanh(A„/i) - i/„(0)]} 



(3.66) 



If we substitute these expressions for the amplitudes A(t)p and B{t)p in the assumed 
form of the particular solution, we obtain 



cos(wt). 



(3.67) 



As a result, we obtain for An{t) the following expression: 



An{t) = An cos(n„t) + Bn sin(n„t) — — Fm cos(a;t) 



(3.68) 



(tj2 - Q?) 

For the amplitude Cm we obtain the same expression as above for An(t), but with 
the index m instead of the index n. Now the potential function can be written in the 
form 



E 



An cos[Unt) + Bn sm(i2„t) - — ^ — — - cos(a;t) 



A Fn \ A Fn 

-—u)cos{u)t)j^Gn{y)j cosh(A„(y + h)) + — w cos(a;t)^ff„(y) sinh(A„(y + h)) 



cos{kznz) 



m=— oo 

A 



Cm cos(n„t) + Dm sin(n„t) - cos(a;t) 



Fn \ A Fm 

^^u) cos{u)t)Y^Gm{y) j cosh(A„(y + h)) + —u)cos{u)t)j^Hm{y) sinh(A„(y + h)) 

(3.69) 



s'm{k^mz) 



which is a function of the unknown constants An,Bn,Cm and Dm- To obtain these 
constants we use the initial conditions for (l){k,y,z,t) given by equations (3.42) and 
(3.43). We obtain 
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= - 62 ^;^^"(^) ~ 62 ^\~^"(^) tanh(A„/i), (3.70) 

Bn =0, (3.71) 

^™ = — 2 _ o2 ^ ?;2 ^\ — - — Hm{0) tanh(A„/i), (3.72) 

Dm =0. (3.73) 
The final form of the potential function (/)(A;, y, z, t) is given by the equation 



cos(n„t) - cos(a;t) cosh(A„/i) . . 

+ u) — (cos(i2„t) — cos(a;t)) 



i 62 L A2(a;2-n2) cosh(A„/i) 



A2 



cj sinh^(A„/i) 
A2 cosh(A„/i) 



cos(n„t) 



, , A sinh^(A„(y + /i)) ^ x ■ /, x 
cosh(A„(y + h)) + -^^Fn cos(a;t) \ cos{k^nZ) 



+ E 



cos(fim^) — cos(wt) cosh(ATO/i) 



1 62 [ A^(a;2-n^) cosh(A™/i) 



+ 0;- 



A2 



(cos(Qm^) — cos(a;t)) 



uj sinh^(ATO/i) 
A^ cosh(A„/i) 



cos(fiTO^) 



A sinh^ (Afniy + h)) , A 
cosh(A„(y + /i)) + T^ioFm -T^ cos(wt) ^ sm{k^rnZ). 

(3.74) 



We are interested in the displacement of the free-surface r]{k,z,t), which is given 
in terms of the potential function (/){k,y,z,t) according to the equation (3.31). Then 
the cosine transform of the free-surface displacement is given in terms of the Fourier 
transform of the potential according to the equation 



1 d(f) 

ri[k,z,t) = — — [x,Q,z,t). 
9 dt 



(3.75) 



If we apply this equation to the expression for (/){k,y,z,t) given by equation (3.74), 
we obtain 
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fj{k,z,t) = 

n=—o 
oo 

+ E 



Ujfln 



(w sin(a;t) — n„ sin(n„t)) 



cos(A;2„2;) 



(cj sin(a;t) — 0^ sin(nmt)) 



cos(A;^„2;)| . 

(3.76) 



3.2 Fourier Integral Solution. 



Here we apply the inverse cosine transform to the expression above for the cosine trans- 
form of the free-surface displacement. The inverse cosine transform is given by equation 
(3.36), and we apply it to the equation (3.76) to obtain the free-surface displacement 



+ 



E 

«=— oo 

oo 



1 r- 



gb^ 



{uj sin(6jt) — VLn sin(0„t)) cos{kx)dk 



(u) sin(a;t) — ilm sin(Om^)) cos{kx)dk 



Sm{k^rnZ) 



(3.77) 



The integrands in the integrals above apparently have poles in the complex k plane 
for wave numbers solutions of co"^ — ^^i^) — ^- ^n{k) approaches icj, we have that 
fln{k) sm{u}t) approaches a;sin(a;t) in the same fashion, so there is no singularity in 
the integrand and the integral is well behaved. To obtain the free-surface displacement 
we evaluated numerically the inverse cosine transforms appearing in equation (3.77). 
Results from these simulations were used to generate animations of the evolution of the 
free-surface displacement due to the action of the wave maker over the fluid. These 
animations are discussed in the next section. 



3.3 Numerical Results. 

Here we show results from the numerical evaluation of the inverse cosine transforms ap- 
pearing in the equation (3.77) for the free-surface displacement. We display the evolution 
of the free-surface displacement in time through the numerical evaluation of equation 
(3.77). We generated animations for the evolution of the free-surface displacement due 
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to the action of the wave maker at a; = 0. Here we discuss the examples and we give 
links for the movies associated with these examples. 

• We consider that the displacement of the wave maker coincides with the first 
cosine wave mode in the z direction. The excitation frequency is above the cut- 
off frequency for the first cosine wave mode. With this type of excitation, the 
only wave mode taking part in the solution is the first cosine wave mode. Since 
the wave maker starts from rest to the harmonic motion, it excites initially all 
wave frequencies and generates a transient which propagates along the channel 
and is followed by a nono-chromatic wave train (the cosine wave mode) with 
frequency equals to the excitation frequency. The transient has a wave front which 
propagates with the maximum group velocity possible for this cosine wave mode. 
For the depth h = 0.1 meters, figure 5 illustrates the maximum group velocity for 
the cosine wave modes. The maximum group velocity possible Cg^max is the group 
velocity of the cosine wave mode with k^n = ^ (n = 0). Then, for a given time 
instant t, there is no wave disturbance at positions x > Cg^maxt- The transient 
for a given instant t stays in the region Cg^maxt > x > Cg{uj)t, where Cg{uj) is the 
group velocity of the excited cosine wave mode at the excitation frequency w. To 
see the animation associated with this example, click here. 

• We consider that the displacement of the wave maker coincides with the second 
cosine wave mode in the z direction. The excitation frequency is above the cut-off 
frequency for the first cosine mode but below the cut-off frequency for the second 
cosine mode. Again, the wave maker starts from rest to the harmonic motion. 
All wave frequencies are excited initially and a transient develops. The transient 
propagates along the channel, and behind it we are left with only the second cosine 
wave mode, which decays exponentially as we go away from the wave maker, since 
at this excitation frequency the second cosine wave mode is evanescent. Again, 
the transient has a wave front which propagates with the maximum group velocity 
possible for the second cosine wave mode. To see the animation associated with 
this example, click here. 

• We consider that the displacement of the wave maker coincides with the first 



22 




wavenumber (radians per meter) 



Figure 5: Group velocity as a function of the wavenumber k for various values of the 
eigenvalue k^n and water depth h = 0.1 meters. The maximum group velocity for the 
first cosine wave mode {Cg^max) is indicated in the figure. Maximum group velocity for 
the second cosine mode also indicated in the figure. 
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cosine mode in the z direction. The excitation frequency is exactly at the cut-off 
frequency. Again, the wave maker starts from rest to the harmonic motion, and 
initially all wave frequencies are excited. A transient develops and propagates 
along the channel. The transient has a wave front which propagates with the 
maximum group velocity possible Cg^max for the first cosine wave mode. Behind 
the transient we are left with the first cosine wave mode, since it is the only wave 
mode excited by the wave maker. The group velocity of this wave mode at its 
cut-off frequency is zero, so there is no energy propagation along the channel after 
the transient part of the solution is already far from the wave maker. Since the 
energy cannot be radiated away from the wave maker, we see the wave amplitude 
growing with time close to the wave maker. The cosine wave mode resonates with 
the wave maker in this case. To see the animation associated with this example, 
click here. 

• Now the wave maker is a liner function in the z direction {F{z) = z). We show the 
evolution of the disturbance due to the action of the wave maker. We consider all 
modes that take part in the solution. We actually consider only a finite number 
of sine and cosine wave modes. As the wavenumber kzm or k^n associated with 
a wave mode increases, its amplitude decreases, so only a finite number of wave 
modes are significant. Again, the wave maker starts from rest to the harmonic 
motion. We have initially a transient which propagates along the channel. It has a 
wave front which propagates with the maximum possible group velocity, which is 
the maximum group velocity for the first sine wave mode. Ahead of the wave front 
{x > Cg^maxi for a given instant t, where Cg^^nax is the maximum group velocity 
for the first sine wave mode) we have no waves disturbance. For a given instant t, 
the transient stays in the region Cg^maxt > x > Cg{uj)t, where Cg{uj) is the group 
velocity of the first sine wave mode at the excitation frequency u. Behind this 
region we have the steady state solution. To see the animation associated with 
this example, click here. 



